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Abstract. Consider the ensemble of real symmetric Toeplitz matrices whose 
entries are i.i.d random variables chosen from a fixed probability distribution 
p of mean 0, variance 1 and finite higher moments. Previous work |BDJ1 IHM] 
showed that the limiting spectral measures (the density of normalized eigenval- 
ues) converge in probability and almost surely to a universal distribution almost 
that of the Gaussian, independent of p. The deficit from the Gaussian distri- 
bution is due to obstructions to solutions of Diophantine equations and can be 
removed (see |MMS| ) by making the first row palindromic. In this paper, we 
study the case where there is more than one palindrome in the first row of a real 
symmetric Toeplitz matrix. Using the method of moments and an analysis of 
the resulting Diophantine equations, we show that the moments of this ensem- 
ble converge to an universal distribution with a fatter tail than any previously 
seen limiting spectral measure. 
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1. Introduction 



1.1. Background. Since its inception, Random Matrix Theory has been a pow- 
erful tool in modeling highly complicated systems, with applications in statis- 
tics [Wis] ■ nuclear physics |Wigl Wig2 Wig3 Wig4 Wig5] and number theory 
|KSlt IKS2t IKeSnj : see [FMj for a history of the development of some of these 
connections. An interesting problem in Random Matrix Theory is to study sub- 
ensembles of real symmetric matrices by introducing additional structure. One 
of those sub-ensembles is the family of real symmetric Toeplitz matrices; these 
matrices are constant along the diagonals: 





( bo 


bi 


b2 ■ 


■ bN-i 


\ 




b-i 


bo 


bi ■ 






An = 


b-2 


b-i 


bo ■ 








\ bi^N 


b2-N 


bs-N ■ 


■ bo 


) 



Initially numerical investigations suggested that the density of the normalized 
eigenvalues was given by the standard normal; however, Bose, Chatterjee, Gan- 
gopadhyay |BCGj . Bryc, Dembo and Jiang |BDJ] and Hammond and Miller |HMj 
showed that this is not the case; in particular, the fourth moment is 2 2/3 and not 
3. The analysis in |HM] shows that although the moments grow more slowly than 
the Gaussian's, they grow sufficiently fast to determine a universal distribution 
with unbounded support. The deficit from the standard Gaussian's moments is 
due to obstructions to Diophantine equations. 

In [MMSj . Massey, Miller and Sinsheimer found that, by imposing additional 
structure on the Toeplitz matrices by making the first row a palindrome, the 
Diophantine obstructions vanish and the limiting spectral measure converges in 
probability and almost surely to the standard Gaussian. A fascinating question to 
ask is how the behavior of the normalized eigenvalues changes if we impose other 
constraints. Basak and Bose [BBlj . Kargin |Karj and Liu and Wang |LWj obtain 
results for ensembles of Toeplitz (and other) matrices that are also band matrices, 
with the results depending on the relative size of the band length to the dimension 
of the matrices. Another direction is that of Basak and Bose |BB2j . where each 
entry is scaled by the square root of the number of times that entry appears in the 
matrix. These ensembles are special cases of patterned matrices governed by a link 
function; see also |BanBot IBHSj . In this paper we explore another generalization 
by studying the effect of increasing the palindromicity on the distribution of the 
eigenvalues. We begin by listing our notation below, and then stating our results 

inca 
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1.2. Notation. 



Definition 1.1. For fixed n, we consider N x real symmetric Toeplitz matrices 
in which the first row is 2" copies of a palindrome. We always assume N to he a 
multiple of 2" so that each element occurs exactly 2^^^ times in the first row. For 
instance, a doubly palindromic Toeplitz matrix (henceforth referred to as a DPT 
matrix) is of the form: 
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[1.2) 



We always assume the entries of our matrices are i.i.d.r.v. chosen from some 
distribution p with mean 0, variance 1 and finite higher moments. The entries of 
the matrices are constant along diagonals. Furthermore, entries on two diagonals 
that are N/T^ diagonals apart from each other are also equal. Finally, entries on 
two diagonals symmetric within a palindrome are also equal. 



To succinctly keep track of which elements are equal, we may introduce a link 
function ip '■ {1, • • • , N}"^ — )■ {1, . . . , N} and new parameters be such that aij = 
b^(^ij), where 



j\ mod 2" 
- j I mod 2^ 



if 
if 



J I mod 2" < iV/2"+i 
J I mod 2" > A^/2"+i. 



:i.3) 



Each N X N matrix A^y in this ensemble can be identified with a vector in M^/^" 

by An o (6o(-4iv), 6i(Ajv), • • • , &7v/2"(-4iv)). 

Before stating our results, we recall the various types of convergence. The 
following passage is paraphrased from |MMSj with permission. Fix n and for each 
integer N let fln,N denote the set of N x N real symmetric Toeplitz matrices 
whose first row is 2" copies of a palindrome. We construct a probability space 
i^n,N,J^n,N,'Pn,N) by Setting 



Pn,iv {{An e nn,N ■■ bi{AN) e [ai, f3i] for 2 e {0, ... , iV/2" - 1}}) 

= n / Pixi)dxi, (1.4) 
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where each dxi is Lebesgue measure. To each & ^n,N we attach a spacing mea- 
sure by placing a point mass of size at each normalized eigenvalu^ Xi{AN): 

f^n,ANix)dx = ^ ^ (5 - ^'^^ ^ dx, (1.5) 



where S{x) is the standard Dirac delta function. We call /in,Ajv the normalized 
spectral measure associated to A^. 

Definition 1.2 (Normalized empirical spectral distribution). Let A^ & ^n,N have 
eigenvalues Xn > ■■■ > Ai. The normalized empirical spectral distribution (the 

empirical distribution of normalized eigenvalues ) Fn'^^^ is defined by 



(1.6) 

As Fn^^^ {x) = iJ,n,AM{t)dt, wc scc that Fn'^^^ is the cumulative distri- 
bution function associated to the measure ^n,AN- 

We are interested in the behavior of a typical Fn^^^ as we vary in our 

ensembles Vtn,N as — )■ oo. Our main results are that for each n, Fn^^^ 
converges to the cumulative distribution function of a probability distribution. 

As there is a one-to-one correspondence between N x N real symmetric palin- 
dromic Toeplitz matrices whose first row is n copies of a palindrome and M^/^", 
we may study the more convenient infinite sequences. Thus our outcome space is 
^n,N = {bo, bi,. . .}, and if a; = {uq, Wi, . . . ) G fi„,N then 

/•ft 

PToh{uJi e [ai, (3i]) = / p{xi)dxi. (1.7) 

We denote elements of fin,N by A to emphasize the correspondence with matrices, 
and we set to be the real symmetric palindromic Toeplitz matrix obtained by 
truncating A = {bo, bi, . . .) to (6o, • • • , &Af/2"-i)- We denote the probability space 

by (i^n,N, -^n.N, Pn,N)- 

To each integer m > we define the random variable Xm,n;N on fi„,N by 



/oo 
x^^dF^^'^/^ix); (U 
-oo 



note this is the m}^ moment of the measure Hn,AM- 
We investigate several types of convergence. 

(1) (Almost sure convergence) For each m, Xjn,n;N Xm,n almost surely if 

P„,N {{A G nn,n ■■ Xm,n-A^) ^ ^m,n(^) as A^ ^ CX)}) = 1; (1.9) 



^From the eigenvalue trace lemma (Trace(A^) = ^ - Af {Am)) and the Central Limit Theorem, 
we see that the eigenvalues of are of order ^/N . This is because Trace(A^) = X]fj=i ^fj^ 
and since each aij is drawn from a mean 0, variance 1 distribution, Trace(A^) is of size N"^ . 
This suggests the appropriate scale for normalizing the eigenvalues is to divide each by \/N. 



DISTRIBUTION OF EIGENVALUES OF HIGHLY PALINDROMIC TOEPLITZ MATRICES 5 



(2) (In probability) For each m, Xm,n;N — ^ Xm,n in probability if for all e > 0, 



as — )■ oo for all x at which Fx^^{x) = P(Xm,n(^) < is continuous. 
Alternate notations are to say with probability 1 for almost sure convergence 
and in distribution for weak convergence; both almost sure convergence and con- 
vergence in probability imply weak convergence. For our purposes we take Xm,n 
as the random variable which is identically Mm,n-, the limit of the average m}^ 
moments (i.e., limAr_j.oo Afm,n;Ar), which we show below uniquely determine a prob- 
ability distribution. 

Our main tool to understand the Fn^^^ is the Moment Convergence Theorem 
(see |Ta] for example) ; while the analysis in |MMS] was simplified by the fact that 
the convergence was to the standard normal, similar arguments hold in our case 
as the growth rate of the moments of our limiting distribution implies that the 
moments uniquely determine a probability distribution. 

Theorem 1.3 (Moment Convergence Theorem). Let {F/v(a;)} be a sequence of 
distribution functions such that the moments 



exist for all m. Let {Mm}'^=i be a sequence of moments that uniquely determine 
a probability distribution, and denote the cumulative distribution function by ^ . 
//liniAT^oo M„,^N = Mm then limAr^oo Fn{x) = ^(x). 

Definition 1.4 (Limiting spectral distribution). If as N ^ oo we have Fn^^^ 
converges in some sense (for example, in probability or almost surely) to a distri- 
bution Fn, then we say Fn is the limiting spectral distribution of the ensemble. 

1.3. Results. Our main result concerns the limiting behavior (as a function of 
the palindromicity n) of the fin,AN generic as N oo. We analyze these 
limits using the method of moments. Specifically, for each Aj^ we calculate the 
moments of Hn,Ai^ by using the Eigenvalue Trace Lemma to relate the rrfi^ moment 
to the trace of A^. We show the average m*^ moment tends to the rrfi^ moment of 
a distribution with unbounded support. By analyzing the rate of convergence, we 
obtain results on convergence in probability and almost sure convergence, which 
we recall below. 

Specifically, our main result is the following. 

Theorem 1.5. (Convergence in Probability and Strong Convergence) Let 

n be a fixed positive integer, and for each N a multiple o/2" consider the ensemble 
of real symmetric N x N palindromic Toeplitz matrices whose first row is 2" copies 
of a fixed palindrome (see (11. ip for an example), where the independent entries 
are independent, identically distributed random variables arising from a probability 
distribution p with mean 0, variance 1 and finite higher moments. Then as N ^ oo 



lim ¥n,n{\Xm,n;N{A)-Xm,n{^)\> e) = 0; 

(3) (Weak convergence) For each m, Xm,n;N Xm,n weakly if 

Pn,N(^m,„;7v(^) < x) ^ P(X„,„(A) < x) 



(1.10) 



(1.11) 




(1.12) 
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the measures fin,AN ( defined in (11.51) ) converge weakly to a limiting spectral measure 
with unbounded support. If additionally p is even, then the measures converge 
strongly to the limiting spectral measure. 

As in other related ensembles, it is very difficult to obtain closed form expres- 
sions for the general moments of the limiting spectral measure. We can, however, 
analyze the moments well enough to determine the limiting distribution has un- 
bounded support; in fact, as the following theorem shows, it has fatter tails than 
previously studied ensembles. 

Theorem 1.6. (Fat Tails) Consider the ensemble from Theorem \ 1.51 For any 

fixed n > 2, the moments grow faster than the corresponding moments of the 
standard normal, and thus our density has fatter tails than the standard normal. 
Specifically, if fn denotes our density and denotes the density of the standard 
normal, for all B sufficiently large we have J"^ fn{x)dx > (j){x)dx (i.e., the 
probability of observing a value at least B is larger for our limiting spectral distri- 
bution that the corresponding value for the standard normal). 

Similar to other papers, our analysis is based on analyzing the contribution 
from various matchings. We are able to analyze in complete detail the case when 
we have all adjacent matchings. Based on numerical investigations and some 
theoretical calculations, we believe that all configurations contribute equally; see 
Conjecture 15.91 for a precise statement. If true, this would allow us to sharpen 
Theorem II. 6[ 

Theorem 1.7. Assuming Conjecture \5.9i if M2m,n denotes the 2m*^ moment of 
the limiting spectral measure of our ensemble for a given n, then 

M2n.,n > (2m -1)!!. (1.13) 

m 

The limiting spectral measure thus has unbounded support, and fatter tails than the 
standard normal ( or in fact any of the known limiting spectral measures arising 
from an ensemble where the independent entries are chosen from a density whose 
moment generating function converges in a neighborhood of the origin). 

The rest of the paper is organized as follows. We ffist establish some basic results 
about our ensembles and the associated measures in ^ We then analyze the even 
moments in detail in ^ and prove our convergence claims in §H We give the proof 
on the vanishing Diophantine obstructions for highly palindromic Toeplitz matri- 
ces and discuss the configurations of the different matchings of highly palindromic 
Toeplitz matrices; our conjecture on the contribution from various matchings is 
discussed in detail in the appendices. While it is difficult to isolate the exact value 
of these moments, we are able to analyze these moments well enough to prove our 
convergence claims and to have some understanding of the limiting spectral mea- 
sure. The situation is different for both the fourth moment for any palindromicity, 
and we determine the exact values in §5.2.11 

While many of the arguments in this paper about general properties of the 
moments and convergence are straightforward generalizations of those in |HMt 
IMMSj . the higher degree of palindromicity creates numerous technical difficulties 
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which must be overcome if we want to compute actual values of the moments. In 
particular, the combinatorics becomes significantly harder, as can be seen by the 
length of §3, which is the heart of this paper. 

It is worth remarking on the role numerical investigations played in our analysis. 
These studies were essential in highlighting the key features and illuminating the 
structure and the combinatorics. Unfortunately the rate of convergence is signif- 
icantly slower than in the ensembles in |HMt IMMS] , and it is thus non-trivial to 
extract useful data from these simulations. We end the paper with a brief dis- 
cussion in Appendix |X] on a fast way to test conjectured answers to the difficult 
combinatorial problems. 

The results in this paper arise by increasing the palindromicity of the matrices, 
which leads to highly patterned matrices. A natural question is what happens if 
instead we weaken the structure. In a sequel work, M. Kologlu, G. Kopp and S. 
J. Miller |KKM] do just that. They study the ensemble of real symmetric period 
m-circulant matrices, where each diagonal is now periodic with period m (a real 
symmetric circulant matrix has constant diagonals that wrap around). This gives 
an interpolation from the highly structured symmetric circulant matrices, whose 
limiting eigenvalue density is a Gaussian, to the ensemble of all real symmetric 
matrices, whose limiting eigenvalue density is a semi-circle. The combinatorics 
can be reinterpreted as a counting problem in algebraic topology, and closed form 
expressions are obtained. The limiting spectral measure is the product of a Gauss- 
ian with a polynomial of degree 2m — 2, and rapidly converges to the semi-circle 
as m tends to infinity. 



2. DiOPHANTINE FORMULATION 

In this section we begin our analysis of the moments. We prove some combina- 
torial results which restrict the number of configurations which can contribute a 
main term; we then analyze the potential main terms in the following section. 

Recall that for each matrix An & ^n,N we associate a probability measure by 
placing a point mass of size at each of its normalized eigenvalues Xi{AN): 

fin,ANix)dx := ^ ^ - ^'^^ ^ dx, (2.1) 

where 6{x) is the Dirac delta function. Thus the /c*^ moment of /i„,Ajv(^) is 

/oo 1 ^ 

i=l 

The expected value of the /c*^ moment of the N x N matrices in our ensemble, 
found by averaging over the ensemble with each A^- weighted by (??) and using 



8 



JACKSON, MILLER, AND PHAM 



the Eigenvalue Trace Lemma, is 



1 




l<ii,...,ik<N 



1 



l<h,...,ik<N 



(2.3) 



where from (??) the expectation equals 




(2.4) 



We let Mfc „ be the limit of the average moments; thus 

Mk,n ■= lim Mk,n;N; 



(2.5) 



we will prove later that these limits exist. 

Our goal is to understand the Mk^n, i-e., the limiting behavior of the moments 
in these ensembles. We use Markov's Method of Moments, which we summarize 
below. This is a standard method for proving results in the subject; a nice explicit 
summary of this method begins Section 3 of |BBlj . 

• We first show Mm,n = limN^oo Mm,n;N = limN^ooE[Mm,n;N{AN)] exists 
for m a positive integer, with the Mm,n's satisfying Carleman's condi- 
tion: ^^=1 Mg^*^^™ = oo. As these are the moments of the empirical 
distribution measures, this implies that the M^ ^'s are the moments of a 
distribution. 

• Convergence in probability follows from analyzing the second moment, 
namely showing Var(Mm„n;Ar(^Ar) — Mm,n) tends to zero as — )■ cxd. 

• Almost sure convergence follows from showing the fourth moment tends to 
zero and then applying the Borel-Cantelli lemma. 

We do the convergence calculations in ^ in this and the next few sections we 
determine the limiting behavior of the ensemble averages. 

The odd moments are readily determined, as counting the degrees of freedom 
show the average odd moments vanish in the limit as A^ — ?■ oo. 

Lemma 2.1. All the average odd moments vanish in the limit; i.e. limAr^oo ^2m+i,n;Ar 



Proof. For the (2m + l)*"* moment, we consider IE[&V'(n,j2)^V'fe,i3) " " " ^i>{i2m+i,h)\'^ we 
may write this as IE[&^^ ' ' ' with ri + ■ ■ ■ + = 2m + 1 and the 6^'s distinct. As 
2m + 1 is odd, at least one is raised to an odd power. If any of these occur to 
just the first power, then the expectation is zero as the 6's are drawn from a mean 
zero distribution. Thus at least one of the 6^'s above occurs at least three times, 
and every hi occurs at least twice. The maximum number of distinct 6^'s occurs 
when everything is matched in pairs except for one triple matching. Thus there 
are at most m different fe^'s in our tuple, and the number of tuples is bounded 
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independent of N. We have two degrees of freedom from the first matching of the 
ftf's and one degree of freedom for each other matching]^ for a total of at most 
m + 1 degrees of freedom. Thus the number of indices ii, . . . ,ik G {1,..., N} that 
can contribute to the moment in (I2.3p for a given matrix is 0„(A^™'"^^) (where the 
big-Oh constant may depend on n, as the larger n is the more choices we have for 
diagonals). As we divide by A^™'+3/2 f|2.3l) . the odd moments are On{N~^^^), 
and thus vanish in the limit as — )■ oo. □ 

Corollary 2.2. For fixed n, as N ^ oo there is no contribution to the average 
2m*'' moment from any tuple where the bi 's are not matched in pairs. 

Proof. The corollary follows from a similar analysis as in Lemma 12.11 □ 

From the above corollary, we see that in order to study the eigenvalues of our 
matrices we need to know how many different ways the k = 2m entries (the 

i _^^'s) in our tuples can be matched into k/2 = m pairs. Letting r!! = r(r — 
2) (r — 4) ■ ■ ■ , where the product stops at 1 if r is odd and 2 if r is even, we see 
there are at most (2m — 1)!! ways to match in pairsB Note (2m — 1)!! is the 
2^th jYioment of the standard normal, and has the combinatorial interpretation 
of being the number of ways of matching 2m objects in m pairs where order does 
not matter. For each legitimate matching we obtain a system of m equations, one 
for each pair of entries, for which the number of solutions is the contribution of 
the matching to the 2m**^ moment. 

In order to understand the even moments, we need to know more about the 
permissible matchings, and how many choices of the indices lead to valid config- 
urations. In the original case of the ensemble of real symmetric Toeplitz matrices 
jHM] . the only way any two entries bg could match was for them to lie on the 
same diagonal or on the reflection of that diagonal over the main diagonal. That 
is, they matched if and only if 

\ip-ip+i\ = \ii-ii+i\. (2.6) 

For highly palindromic Toeplitz matrices, more relations give matchings (as seen 
in the investigation of palindromic matrices in |MMS] ). An entry for which the 
absolute value of the difference between its indices is in a given congruence class 
modulo 2" can match with another entry if and only if it is in the same congruence 



For example, say &i/)(ii,i2) = ^^/>(i„,i„+i)j with ii our first index. Both ii and i2 are free 
variables and we have N choices for each; however, iy is not (it will have occurred in a matching 
before this point), and iy+i is determined by requiring the two bis under consideration to be 
equal. The number of choices for i^+i depends on n (the larger n is, the more diagonals work); 
what matters is that the number of choices for i^+i is bounded independently of N. Whenever 
we have a new pair, we have a new choice for the value of the link function, and thus gain a 
degree of freedom. 

"'There are (^™) ways to choose the first pair, (^™2^^) ways to choose the second and so on; 
we must divide by ml as it does not matter which pair we call the first. The claim follows by 
elementary algebra. Alternatively we can prove this by induction. Assume there are {2m — 3)!! 
ways to match 2m — 2 objects in pairs. If we have 2m objects, there are 2m — 1 choices of an 
element to pair with the first element in our list, and then by induction there are (2m — 3)!! ways 
of pairing the remaining 2to — 2 elements. 
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class or its negative. That is, two entries ai^ip^^ and a^^^-^^ can be matched in a 
pair if and only if their indices satisfy one of the following relations: 

(1) there is a Ci G {(-[^^^^J + k-l)§-l\ke {1,...,2"}} such that 

\ip-ip+i\ = \ii - ii+i\ + Ci, (2.7) 

(2) there is a C2 G {(L^^^i^J + k)§ | A; G {1, . . . ,2"}} such that 

lip - ip+il = -Ih - ii+il + C2; (2.8) 
as is standard, [x\ represents the largest integer at most x. 

As a consequence of 02.71) and fl2.8p . for the matchings above there is some C 
such that 

ip - ip+i = ±(«« - ii+i) + C. (2.9) 
As there are two choices for sign for each of the m matchings, there are potentially 
2™ cases that can contribute. We now prune down the number of possibilities 
greatly by showing only one case contributes in the limit, namely the case when 
all the signs are negative. 

In the Toeplitz ensembles studied in |HM] and |MMS] . it was shown that any 
matching with a positive sign (i.e., as in (12. 9p ) in any pair contributes a lower order 
term to the moments, and thus it sufficed to consider the case where only negative 
signs occurred. A similar result holds here, which greatly limits the number of 
cases we need to investigate. Note by Lemma 12.11 we need only investigate the 
even moments. 

Lemma 2.3. Consider the contribution to the 2m^^ moment from all tuples [ii, . . . , 
^2m) in which the corresponding h^'s are matched in pairs. If an ai^i^^-^ is matched 
with an a,^^^^ with a positive sign (which means 

ig - iq+i = +{ii - ii+i) + C 

for some C as defined in (12. 7p or (12. 8p ). then this matching contributes Om{^/N) 
to M2m,n;N (J'lT'd thcrcforc the contribution of all but one of the 2™ choices for the 
m signs vanishes in the limit as N ^ 00, with only the choice of all negative signs 
being able to contribute in the limit. 

Proof. The argument is essentially the same as in |MMSj . Briefly, the idea is that 
if there is ever a positive sign then we lose a degree of freedom, leading to a lower 
order contribution. 

For any tuple {ii, . . . , i2m) in which the corresponding bes are matched in pairs, 
there exist k equations, one for each pairing, of the form 

iq-iq+i = ei{ii - ii+i) + Ci where ei = lor -1. (2.10) 

Let Xi,X2, ■ . . , X2m denote the absolute value of the difference between two indices 
of each entry (so for cii^^i^^^ it would be xi = \ii — ii+i\), and let Xi = ii — 12^X2 = 
12—12,, ■ ■ ■ and X2m. = i2m—ii (i-e., the unsigned differences). It follows immediately 
that 

2m 

^x, = 0. (2.11) 

1=1 
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Each be expressed in two ways. By breaking the absolute value sign in 

(12 .yp or fl2.8p . we have Xr = rjjXj for some j with 77^ = 1 or — 1. We can also 
express it through an equation like the one in fl2.1ip such that Xr = e^xi + Cr for 
some /. Thus 

Xr = r]jXj = erXi + Ci. (2-12) 

Then since = 1, 

Xl = erTjjXj — trCi. (2.13) 

Note each absolute value of a difference occurs twice, as everything is matched in 
pairs. We therefore have 

2m m m 

= "^[qjXj + {erVjXj - erCn)] = '^{njXj{l + er) - erCj) = 0. (2.14) 

i=l j=l j=l 

If any = 1, then the Xj's are not linearly independent and we would have less 
than m + 1 degree of freedomj^ The contribution from such tuples to the moment 
in 02.30 for a given matrix is therefore 0{1/N) (as we divide by A^'"^^), which 
vanishes in the limit as — > 00 and can thus be safely ignored. □ 

Lemma 12.31 immediately implies 

Lemma 2.4. // the indices of ai^ii^^ and ai^i^_^_-^ satisfy (12. 7p for some Ci, then 
\iq - iq+i \ = \ii - ii+i \ + Ci implies 

,-iq+i = -{ii - ii+i) + Ci ^^\ig-ig+i = -{ii - ii+i) - Ci 
iq > max{ig+i,ig+i + Ci} ]^iq < min{ig+i, i^+i - Ci}. 

(2.15) 

Similarly, if the indices satisfy (12.80 for some C2, then \iq — iq+i\ = — + 
C2 implies 

iq-iq+i = -{U - k+i) + C2 or <'^~^'^'^^ ^ -{U-ii+i)-C2 

iq+l < iq < iq+l+C2, ^iq+l - C2 < iq < iq+l- 

(2.16) 

Instead of considering each value of C (either Ci or C2) individually, we will 
consider a pair of constants Ci, C2 such that C1 + C2 = N — 1. We claim that this 
removes some of the Diophantine obstructions that arise when evaluating fl2.15p 
or (12.161) individually. Given an entry we can associate each value of C with 

one diagonal whose entries, generally denoted by ai,,i,.+i, all equal Except for 

the main diagonal, every other diagonal has fewer than entries and therefore 
the index v G {a, . . . , 6} where 1 < a < b < N rather than v € {1, . . . , N}. 
Here we only need to restrict one of the two indices of aj^j^^^ and the other one 
will automatically be determined. However, by considering ai^i^^-^ on a pair of 
diagonals associated with Ci,C2, we can take the index v (or v+i) to be any 
value between 1 and N. Furthermore, except for 0(1) values, the first index of 
entries from the pair of diagonals associated with Ci, C2 are distinct, and similarly 



^As in the proof of Lemma 12.31 the first pair gives us two degrees of freedom and each 
subsequent pair gives at most one degree of freedom. If the Xj s cirG not linearly independent, 
there can be at most to — 1 independent x/s, and thus at most to degrees of freedom. 
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for the second index. Therefore, if aj^j^_^^ is on the diagonal associated with Ci 
and a./^' is on the diagonal associated with C2, then for some a, 6 G {1, . . . , N}, 
we have: 

ve{a,...,6} 

ie{0,...,a}U{b,...,N} ^ ' ' 

3. Properties of the Even Moments 

In Lemma 12.11 we showed that the average odd moments vanish in the limit. 
In this section we analyze the even moments. While the low moments may be 
computed by brute force (we provide the computation for the fourth moment in 
^ and discuss its consequences), similar to other ensembles we are unable to 
obtain nice closed form expressions for the higher moments in general, although a 
combination of numerical simulations and some partial results suggest the answer 
for the doubly palindromic case; see Appendices \M and [Bl 

3.1. General Properties. We first handle the zeroth and second moments, and 
then turn to the higher moments. 

Lemma 3.1. Assume that p has mean 0, variance 1 and finite higher moments, 
and fix the degree of palindromicity n. Notation as above (see (12.21) and (12. 5p ). for 
all An we have Mo,„;Ar(Ajv) = 1 and M2,n;7v(^iv) = 1? which implies the average 
moments in the limit are both 1 (explicitly, Mo,„ = 1 and M2^n = 1/ 



Proof. From (12.21) . we see Mo,n,Af(^7v) = 1- For the second moment, we have 

1 

iV2 



l<ii,i2<N 

l<ii,i2<N l<ii,i2<N 



Since we choose the 6's from a distribution with mean zero and variance 1, the 
expected value above is just the variance (which is 1), and hence M2,n;N = 1, 
which implies M2^n = lim7v-s>oo ^2,n;Af = 1- D 



We now consider the general even moments. By Corollary 12. 2[ the only contri- 
butions to the moments M2m,n;N (see (12. 3p ) that survive as — )■ 00 is when the 
matched in pairs. There are (2m — 1)!! such matchings; we need to 
determine the contribution of each matching to M2m,n;N- 

Each of the (2m — 1)!! matchings, hereafter referred to as a configuration, leads 
to a system of m equations of the form (I2.15P or (I2.16P (with the C's coming from 
(12. 7p and (12. 8p ). for which each distinct solution gives us one possible choice for 
the tuples (ii, . . . , i2m) and contributes one to the sum. The analysis is completed 
by counting how many valid configurations there are (or at least determining the 
main term). 

Determining the exact value is complicated by the fact that there are many 
ways for an ai.i.^^ and an dj^j^+i to be paired; they must correspond to the same 
&£, but there are many diagonals each can lie on (with the number of diagonals 
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growing with n). Fortunately, we can obtain a weak bound depending on n that 
nevertheless suffices to prove the existence of a limiting spectral measure. By 
standard arguments, it suffices to show the average even moments converge as 
iV — !■ oo to a sequence satisfying Carleman's condition, and then perform a similar 
analysis on the variance (for convergence in probability) or the fourth moment (for 
almost sure convergence). We leave the convergence issues to §11 and instead prove 
the existence of the limits. 

Lemma 3.2. For fixed n, M2m,n exists and 

M2m,n = lim M2m,n;N < (2 " 2")"^(2m - 1)!!, (3.2) 

N^oo 

which implies the M2m,n satisfy Carleman's condition. 

Proof. Fix n and m. Consider one of the (2m — 1)!! pairings where the {aj,j,^^}^Z!\ 
are matched in pairs. We have m equations. For our system of equations we must 
choose m values for the C's, and in each equation by (12.71) and (12.81) there are 
at most 2 ■ 2" choices for a C. By our previous analysis, at most m + 1 of the 
indices ii, . . . , i2m are free. We must analyze the contribution from each choice of 
the indices in (12. 3p . Assume ffist each aj^jj^-^ is matched with a unique ajpi^^-^. In 
this case, the contribution of this choice of indices to (12. 3 p is just the product of 
m copies of the expected value of the second moment of the probability density p; 
as the second moment is 1, each of these adds 1 to (12. 3p . and there are clearly at 
most N^^^ choices of indices. 

If all of the aiiii^^ are not uniquely matched with another ajpi^^^, then it is pos- 
sible to have a larger contribution than 1 to (12. 3p . as the product of the expected 
values could involve fourth, sixth, eighth, 2m*'^ moments. Let Pm be the max- 
imum of the absolute values of the first 2m moments of p. The contribution in 
this case is at most pj^; while this is growing with m, the number of indices that 
can contribute this is at most 0{N^) by our earlier arguments, as we showed we 
lose at least one degree of freedom when items are matched in more than pairs. 

As we divide by A^™+^ in (12. 3p . we find 

M2..,n,N < ^ j^, ^ [Ar-+i . 1 + O(Ar-) ■ pZ] . (3.3) 

Taking the limit as — ?■ oo yields 

M2„,„ = lim M2^,„.;v < (2-2")™-(2m-l)!!. (3.4) 

The existence of the hmit is proved analogously to |BBlt IBDJt IHMt IMMS] (see 
for example Theorem 2.6 of |HM] ) : now that we know M2m,n is bounded, it is easy 
to see that the main term of the contribution from each possible configuration is 
independent of A^. 

It remains to show that the M2m.,n satisfy Carleman's condition by showing the 
sum of the reciprocals of their 2m*'^ roots diverge. Trivial estimation suffices. As 
(2m - 1)!! < (2m)2"^, we have (2m - l)!!-^/^™ > i/2m, and thus 

m m 
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The latter sum is a multiple of the harmonic sum and diverges, completing the 
proof. □ 

4. Convergence 

In §2] and §31 we showed that the limit of the average moments exist as — )■ oo, 
culminating in Lemma where we proved that the moments grow slowly enough 
to uniquely determine a probability distribution. We now show convergence in 
probability, and if p{x) is even we prove almost sure convergence. As these ar- 
guments closely follow those in [HMt IMMS] , we concentrate on the novelties in- 
troduced by the higher palindromicity. We conclude by obtaining lower bounds 
for the moments. These bounds imply that our limiting distributions have un- 
bounded support and fatter tails than the standard normal (possibly the fattest 
tails observed from a random matrix ensemble arising from entries chosen inde- 
pendently from a distribution whose moment generating function converges in a 
neighborhood of the origin). 



4.1. Convergence in probability. We will prove our probability measures con- 
verge in probability. We follow the arguments in [HMt IMMS] . We begin by 
defining our random variables. Let A be a sequence of real numbers to which we 
associate an N x N real symmetric Toeplitz matrix with 2" palindromes, which 
we denote by A^. Thus we may view A as (6o, &2, • • • )' we form An by 
considering the initial segment of length iV/2"+^, taking that as the first half of 
our palindrome, and then building the matrix by having 2" palindromes in the 
first row. 

Let Xm,n;N{.A) be a random variable that equals the m^^ moment of A^ (so 
Xm,n;N{,A) = Mm,n;N{A]\f)) , and Set Mm,n;N to the w}^ moment averaged over the 
ensemble as above (so Mm,n;N = '^[Xm,n;N])- 

Thus, we have convergence in probability if for all e > 

lim Pn({A G ■■ \Xm,n-N - > e}) = 0. (4.1) 

Using Chebyshev's inequality and the fact that Var(r) = E[(r-E[r])2] = E[Y^]- 
E[F]2, we have 

Pn({A G i^N : \Xm,n;N ~ IE[Xm^„.jv]| > c}) 
IE[(Xm,n;Ar — Afm,n;Af)^] 



< 



e2 



N 



£2 



(4.2) 



Thus, it suffices to show 



hm (E[X^^„^^] - Ml^.^^) = (4.3) 



to prove convergence in probability. 
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We have 



l<ii,...,im,<N 

l<jl,-,im<Af 



l<ii,...,im.<Af 



X 5Z E[%i-J2|---%,.-ii|]- (4.4) 

l<jl,-Jm<Af 

We can break this up into two cases. If the entries of the i diagonals are entirely 
distinct from those of the j diagonals, then the contribution to E[X^ „.^] and to 
^m,n;N '^^^^ clearly be the same. Thus, we need to approximate the contribution 
from the cases where there are one or more shared diagonals. The degree of 
freedom arguments of |HM] immediately apply here, though our big-Oh constants 
will now depend on the value of 2" as we now have many more C-vectors to which 
we apply these arguments. Thus, as — )■ oo these two quantities will converge, 
and convergence in probability (and thus also weak convergence) follow. 

4.2. Almost Sure Convergence. We assume that p{x) is even for convenience, 
and use the same notation as above; in particular, 

M™,„ = lim M,n,n;N = lim E[X„,„.;v(^)]. (4.5) 
Almost sure convergence follows from showing that as — )■ oo the event 
{A G 11n : lim A„„.Ar(y4) M„„} 

N^oo 

occurs with probability one for all non-negative integers m. 
By the triangle inequality, we have that 

We have already shown that limAr^oo \Mjn,n;N — ^m,n| = 0, so we need only show 
that \Xm,n]N{A) — Mm,n;N\ almost surely tends to zero. Clearly E[Xm,n;Af(^) — 
Mm,n;N] = 0, and we can modify the arguments in |HMj to show that the fourth 
moment of Xm,n;N{A) — Mm,n;N IS 0^,2" (]^)- All of the degree of freedom argu- 
ments can be applied directly for each C-vector. 

However, Theorems 6.15 and 6.16 of |HM] require greater care as these use more 
than degree of freedom arguments. Fortunately, equations (50) and (51) in |HM] 
hold for any of our C-vectors, so a similar result holds in this case. We then apply 
Chebyshev's inequality to find 

~i - Ar% 

Finally, applying the Borel-Cantelli Lemma shows that we have convergence 
everywhere except for a set of zero probability, thus proving almost sure conver- 
gence. 
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5. Adjacent Matchings, Even Moments and the Tail 

As any distribution with finite mean and variance can be normalized to have 
mean and variance 1, if the distribution is even then the fourth moment is 
the first moment to show the 'shape' of the distribution, and thus merits special 
consideration. We analyze the fourth moment in detail below. We first prove that 
the adjacent and non- adjacent matching configurations contribute equally. We 
then compute the contribution from the adjacent case in §5.2.11 We are able to 
compute the contribution from the adjacent case in the doubly palindromic case 
for any moment; if all configurations contributed equally (which we believe is the 
case) then the 2m*^ moment would just be (2m — 1)!! times the contribution from 
the adjacent configuration. 

5.1. Fourth Moment Configurations. 

Lemma 5.1 (Equal Contribution - Fourth Moment). The non-adjacent configu- 
ration and the adjacent configuration contribute equally to the fourth moment. 




Figure 1. The adjacent and the non-adjacent configurations of the 
fourth moment. By relabeling we see the first and third are equiva- 
lent. 



Proof. For brevity and notational simplicity, we will only present the calculation 
for the non-adjacent case. The adjacent case can be analyzed analogously and 
yields the same main term. That is, we consider the configuration with the fol- 
lowing matchings: 

O'ij = O'kl ^^ 

From the above system of equations relating the matchings, we obtain the corre- 
sponding system of equations for the indices: 

±1^-^'+^ (5.2) 
±\l-t\+B. ^ ' 

Applying Lemma 12. 3[ we need only consider the case where 

J = -(k-l)+A' 

.k = -{l-^) + B' ^^-^^ 
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where either A' = A or A' = —A, and similarly for B\ depending on how the 
absolute value equations resolve. Moreover, we see that 



A + B 



i — j + k — l+j~k + l~i 



0. 



(5.4) 



This implies that A and B must be of the same form, either as in 02.71) or f l2.8p . If 
A is of the form C2 in f l2.8p . then it follows immediately that A = B, whereas if A 
is of the form Ci in (12. 7p . then it can either be that A = B or A = —B. For each 
A we have a system of two equations with four unknowns so we can always pick 
at least two free indices. For convenience, we specify i and j as these free indices 
by choosing aij. Moreover, we assume that we only pick Ojj in the lower diagonal 
half of the matrix so that i > j. By the symmetry of the matrix, picking a^j in 
the upper diagonal half would follow contribute equally. 

We first consider the case where A is of the form C2 in (12. 8 p , and thus A = B = 
C2 for some C2. We then find 



— U 



k\ 



3\+C2 

-l\+C2 



(k-l 

k > I 
J - k 



Co 



(5.5) 



i - Co < I < i. 



We now consider A of the form Ci in (12. 7p where Ci + C2 = N — 1. The value 
Ci is unique for each choice of C2 and the contribution from the pair (Ci,C2) 
complements nicely one another as we show below. We see that A = ±B = Ci. 
We then have 




\t-j\ + Ci 
\l-i\±Ci 



(k-l = 

k < I 

j-k = -{l-t) + 
I < i or I > i + Ci. 



-ii-j)-Ci 



(5.6) 



Since we have already picked the first entry aij, we are left to choose the entry 
ttki- Our choice of Ci (or complementary C2 since the pair is unique) indicates 
the diagonals that an lies on, which gives us the restrictions on i. Finally, as 
only one of the indices k or / need to be specified (since the other is restricted by 
the diagonal), without loss of generality we choose /. We now use our previous 
analysis from (I2.13P and Lemma [23] to analyze the diagonals associated to A = Ci 
and A = C2- Since except for the main diagonal, every other diagonal has less 
than entries, if we choose aki on the diagonal A = C2, then there exist integers 
a,b E {!,..., A^} such that a < I < b and either a or 6 must be A^. If we choose 
aki on the diagonal A = Ci, then I < a or b < I. 

(1) On the diagonal associated with A = C2: 

I e {i-C2,...,i}n{a,...,b}. (5.7) 

(2) On the diagonal associated with A = Ci. 

/e ({0,...,^}U{^ + Ci,...,Ar})n({0,...,a}U{6,...,Ar}). (5.8) 
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Therefore, there are exactly C2 out of + 1 values of i we can pick (or exactly 
Ci out of + 1 value of i we cannot pick). Since we have A^^ choices for picking 
the initial entry ajj, the contribution to the fourth moment from the pair (Ci, C2) 
is given by 



+ k 



(5.9) 



This contribution only depends on the initial choice of aij and the choice of A. 

Repeating this analysis for the adjacent case and summing over all possible 
choices of A of the form Ci, we obtain the same contribution to the fourth moment 
from either configuration. □ 



5.2. The Adjacent Case. In this section we analyze the Adjacent Case in de- 
tail, as this configuration is easier to study and more easily generalized. Since the 
contributions to the fourth moment from the adjacent and non-adjacent configura- 
tions are identical, this allows us to calculate the fourth moment for any number of 
palindromes. We can also calculate the contributions in the general Adjacent Case 
for any even moments for the doubly palindromic Toeplitz matrix. In principle we 
could use the same ideas to calculate by brute force the adjacent configurations of 
any even moment of an ensemble with a greater number of palindromes, but we 
could not find a closed form expression for these moments. We conjecture that all 
configurations contribute the same main term in the limit, and provide numerical 
support in Appendix lAl 

5.2.1. Determining the Fourth Moment. For the fourth moment, we have four 
indices i,j,k, and /, and we consider an adjacent matching where 

Cljj djk: ^kl ^li- 

We think about this as follows. A pair i and j gives us a matrix element Ojj; we 
want to find all pairs j and k such that aij = ajk- This could happen by having 
the two on the same diagonal, or it could happen that ajk is on a palindromically 
equivalent diagonal. As the formula for the fourth moment of our matrix An 
involves division by N^, we need only worry about situations where we have on 
the order of tuples. Clearly we may choose i and j freely. The matching then 
forces there to be on the order of 1 choice for k (the exact answer depends on n, 
the degree of palindromicity; what matters is that the answer is independent of 
N in the limit), and on the order of N choices for I. The last is important, as 
unless the number of choices of / is proportional to N, we will obtain a negligible 
contribution from the matching aij = ajk and aki = an. Exploiting the symmetry 
of the matrix, this reduces to choosing k so that aij = akj and a^ = an. That is, 
in addition to matching aij and akjwe want row i and row k to match well. 

We isolate some of the most useful features of our matrices in the following 
lemma. The proof follows immediately from the previous discussions and the 
structure of the matrices in our ensemble. 



Lemma 5.2. Fix n and consider the ensemble of N x N real symmetric palin- 
dromic Toeplitz matrices with 2" palindromes in the first row. The main diagonal 
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Figure 2. An example highlighting matchings for / in medium 
shading, with mismatching in dark shading. Note that any anoma- 
lous matchings won't contribute in the limit. 

is the only place (excluding the border of the matrix) where bo occurs once rather 
than twice. This implies the following useful properties. 

• Moving to the corresponding point in the next palindrome can require ei- 
ther moving N2~"' — 1 elements when crossing the main diagonal or N2~"' 
elements otherwise. 

• Let c G {1, 2, ... , 2^" — 1}- pictured in Figure 0, a given row and the 
row cN2^"' rows down from that given row do not match perfectly, but 
rather become unaligned when one row has reached the main diagonal but 
the other row has not. Moreover, the row cN2~^ — 1 rows down starts out 
unaligned, but then becomes aligned in this same region. Furthermore, only 
rows of this form match up well with the original row. 

Proof. The first item follows directly from the observation that the main diagonal 
is the only place where bo appears once rather than consecutively. We also see that, 
neglecting the first row which starts on the main diagonal, that the first elements 
of a row and one cN2~'^ rows away match initially. Moreover, they evolve the 
same way when moving from left to right, except when the first one hits the main 
diagonal, in which case it skips forward one place in the palindrome, in which case 
they do not match except possibly for repeating elements at the beginning/end 
or middle of palindromes, like bo- However, once the second row hits the main 
diagonal, it also skips forward, and they become realigned. The case for rows 
cN2~^ — 1 rows away from each other is argued similarly. 

To prove that no other rows match sufficiently well we need to show that there 
are only 0„(1) matchings in any of the other rows. Suppose we do have a matching 
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in one of the other rows. Since we can't be at the corresponding point in the 
pahndrome, we must be at the other end of the pahndrome. Unless these are 
the special repeating elements at the beginning or middle of a palindrome these 
two rows will evolve differently, so although there may be additional anomalous 
matchups, there will certainly not be more than four per palindrome, giving us the 
desired maximum of 0„(1) possible matchings. If they are the special repeating 
elements, then the rows can match up well, but in this case there are at most five 
such elements per row, so we again have a lower order term. □ 

With this lemma in hand, we can now calculate the contribution from a specific 
constant, which will then allow us to calculate the contribution from the adjacent 
configuration. 



Lemma 5.3. Let c E {0,1, ... ,2"' ^} and k = i + cN2 There are then 

N^ + On{N^) (5.10) 



2" - ^ 



good matchings, whereas if k = i + cN2 " — 1, then there are 

On„2 „3 

^^5^iV3 + 0„(iV2) (5.11) 

good matchings, where the big- Oh constants depend on n (which is fixed)^ 

Proof. We begin by noting that by Lemma 15.21 above, choosing k so that Oij and 
Ofcj are at corresponding points in a palindrome guarantees that = a^j and 
that there are On{N) choices of I satisfying aki = an, as desired. Moreover, if 
aij and Okj aren't at corresponding points in the palindrome, then there are only 
0„(1) good choices of /, and since there are at most n such possible cases, this 
contribution can be ignored. Thus, we only consider the cases where Oij and akj 
are at corresponding places in a palindrome. 

We now consider the case when k = i + cN2~^, hence a^j and a^j must be 
on the same side of the main diagonal in order to match. Moreover, to have 
k G {1,2,..., N} we must have i G {1, 2, . . . , N — cA^2~"}. Another restriction 
arises from the fact that they are on the same side of the main diagonal. We note 
that we won't cross the main diagonal when moving down from any Ojj below 
the main diagonal to akj. There will similarly be no crossing if aij lies more than 
cN2~^ elements above the main diagonal. This defines two right-triangular regions 
of height N — cN2~^ + 0„(1), which in total gives a square of area 

2 

N'^ + On{N) (5.12) 



from which to choose Oij, thus giving that many valid choices of a^j. We also have 
the restriction on the values of / as explained in Lemma [5^ leaving — cN2~'"' + 



^The constants may be taken to depend on c as well; however, as n is fixed and c G 
{0, . . . , 2"~^}, we may take the maximum of all the constants and may replace c dependence 
with n dependence. 
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0„(1) good values of / for each of these aij. In total cN2 " contributes 



c 



N-^ + On{N') (5.13) 

matchings to the fourth moment. 

Next we consider the case where we cross the main diagonal when moving from 
ttij to akj, so that k = i + cN2~^ — 1. In this case, the area of values of aij from 
which we will cross the diagonal to give a matching will be mostly defined by 
the parallelogram bordered by the triangles from the previous constant. However, 
there may also be additional strips as depicted in light shading in Figure [3l but 
these will only be of width 1, so the area is essentially that of the parallelogram 
of height - cN2-'' + 0„(1) and width cA^2-" + 0„(1). There will also be 
cA^2~" + 0„(1) good values of Z, so in all this constant contributes 

^ ^ N^ + On{N^) (5.14) 
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Figure 3. Regions where k 
in medium shading, whereas those where k 
in dark shading. Regions where both are satisfied are indicated in 
light shading: These are 1— dimensional, and thus won't contribute 
in the limit. 



We argue similarly for the negative constants {— cA^2^", — (cA^2~" — 1)} for 
cG {1,2,...,2"' — 1}, in which case we are now moving up c palindromes, and 
either crossing or not crossing the mian diagonal, respectively. We easily see that 
this is essentially switching the roles of a^j and aj^, so the contributions will be 
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the same. If we repeat the analysis above we find regions of identical size that 
thus give identical contributions to the fourth moment. Pictorially, what happens 
for a negative constant is that of the positive one rotated 180 degrees. Thus, the 
contribution to the fourth moment will be given by the contributions from the 
positive constants ({cA^2-", cA^2-" - 1} for 6 G {1, 2, . . . , 2" - 1} multiplied by a 
factor of 2 to account for the negative constants. 

Theorem 5.4. Fix n and consider the limit as N ^ oo of the average fourth 
moment of our ensemble. The contribution from one of the adjacent matching 
configurations (i.e., aij = ajk and = an) to this limit is 

^2" + l2-". (5.15) 
3 3 ^ ' 

Since all configurations contribute the same main term by Lemma \5. 1\ we have 

M^^n = 2"+i + 2-", (5.16) 
which is asymptotic to 2^^^ as n — )■ cxd. 

Proof. For each value of c, we note that the contribution to M4^„(A^) is 

2"-c\^_, /2"c2-c3 
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2" - c\ ^ 2"c2 - ^ / 1 



+ ^^^ + On{j^]. (5.17) 



Taking the limit as — i- oo yields the contribution to M4 „ is 



2"-c\^ 2"c2-c3 



+ 5 . (5.18) 

We sum over all values of c, multiply by 2 to account for the negative constants, 
and include the contribution from C = 0, known to be 1 from |MMS] to obtain 
the contribution to the fourth moment from the adjacent matching case: 

= l + ((2" - c)' + 2V - c') . (5.19) 

c=l 



Extending the sum above to include c = cancels the first and last terms of the 
sum, but we must subtract 4 to compensate. This then leaves a sum of squares 
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which is easily evaluated: 



2 



23n 

c=0 



-1 + 
-1 + 



2 2"(2" + l)(2-2" + l) 

22^ 6 

(l + 2-")(2-2" + l) 



3 

= -1 + ^(2 -2" + 2 + 1 + 2-") 

= -2" + i2-". (5.20) 
3 3 ^ ' 

Multiplying by three to account for the two adjacent configurations and one non- 
adjacent configuration then yields 

M4,„ = 2"+i + 2~", (5.21) 

the desired result. □ 

5.2.2. Adjacent Matching Case for Doubly Palindromic Toeplitz Matrices. To cal- 
culate the contribution from the case of all adjacent matchings to even moments of 
doubly palindromic Toeplitz matrices, we simply generalize the pictorial method 
from above. For the 2m^^ moment, we find that our final system of equations for 
the adjacent configuration becomes 

is = ii + Ci 
H = H + C2 = ii + Ci + C2 



H — i2m-l + Cm — il + Ck- 

k=l 

Remark 5.5. The even indices don't appear because the n^^ matching gives the 
equation i2n — i2n-i = —{i2n+i— i2n) + Cn, and the i2n terms cancel. However, for 
each non-zero constant Ci, we will have a picture similar to FigurelE, which limits 
the number of good values of the even indices i2i analogous to the restrictions 
on I for the fourth moment. Moreover, as each i2n+i is related back to ii, the 
difference between the maximum and minimum partial sums must be strictly less 
than N + 0(1) or we lose a degree of freedom. 

Before deducing the main term in the case of all adjacent matchings, we first 
set our notation. 

Definition 5.6. A C-vector is the ordered collection of constants relating the 
odd indices to each other. In the example at the beginning of this subsection, the 
C-vector would be (Ci, C2, . . . , Cm)- A core of a C-vector is the ordered collection 
of nonzero constants in the C-vector. That is, we collapse down the C-vector to 
its core by removing all of the zero constants from it. We can also think of the 
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C-vector as being built up from the core by adding back the zeros in the correct 
places. 

Theorem 5.7. The contribution of one adjacent configuration to the 2m}^ moment 
averaged over the ensemble of doubly palindromic Toeplitz matrices equals 

-2 + 2-™ 6™ (5.22) 

Proof. The foUowing observations greatly simplify the analysis for this case: 

• If the constants iy appear in the C-vector (Ci,. . . ,C„), then iy — 1 
cannot occur as a main term. If it did, we would lose a degree of freedom 
in 12, as aj^j2 would need to lie on a very specific set of diagonals. 

• If some Cj is non-zero, then the next non-zero C chosen must be —Cj, as 
we would otherwise lose a degree of freedom in ii. 

Now consider the 2m*'^ moment, which will have a C-vector of length m. We can 
then consider a subset of length k (k even) of (y,— y,...,^,— y) that forms the 
core of the C-vector, with the remaining entries being zero. There are then ("^j 
distinct ways to insert the remaining zeros, and thus ('^) ways to build a C-vector 
around this core. 

We now consider the contribution from each of these C-vectors. By Remark 

N 
2 



\5.5\ there will be y values of ii to choose from, and there will be k other i2i 



(corresponding to the k nonzero C/) that will have (A^ — y) + 0(1) good values. 
Thus the contribution for each of these cases will be (l)'^^^- The total contribution 
to the 2m*^ moment from this configuration, summing over all possible C-vectors, 
is therefore 

fc+i 

(5.23) 



E 

k even 
k=2 



m\ / 1 

k U 



If we pull out a factor of | and include m = in the sum, we can use the binomial 
theorem to express this as 

^^^,,l]\(,.lY).l (5.4) 



4VV 2 V 2 2 



The contribution from a core of (— y,y,...,— y,y) will be the same as that 
above. The cores of (±(y — 1), =F(y — 1), • • • , ~ -'-)' "F(t ~ ^)) ^le similarly 
analyzed, and they will also have the same contribution since N — ^ + 0{1) = 
Y -|- 0(1), so we multiply fl5.24p by 4. We also include the contribution from 
the 0- vector, which is 1 for the adjacent case. Thus, the contribution from each 
configuration is 

1 + + 1- + (^1 _ = -2 + 2- 6- j , (5.25) 
completing the proof. □ 
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Remark 5.8. Unfortunately, this method does not readily generalize to matrices 
with a greater number of palindromes. The fundamental reason is that the obser- 
vations made at the beginning of Theorem \5. 7| no longer hold for these matrices, 
which then makes it tremendously more difficult to generate all valid C -vectors. 

To demonstrate these difficulties, we investigated the G*'' moment of a matrix 
with four palindromes. While we can construct C -vectors in the same manner 
as for the doubly palindromic ensemble, we clearly will be missing a substantial 
number of possible C -vectors. For instance, we miss the vector (y, ^, — 

In addition to these vectors, we have even more problematic vectors such as 
(y' T ~ ~ ~ ■^))' ''^hich turns out to be valid for chosen within a certain 
parallelogram shaped band of the matrix. These new vectors, in which '^mixing" 
is important, are hard to systematically account for, making it quite difficult to 
determine precisely which C -vectors to include and which to exclude. While we 
could in principle calculate the adjacent contribution to any given moment for any 
number of palindromes, there is no apparent method that will simultaneously work 
for all of these possibilities. Similar to other investigations on related ensembles, 
we are left with general existence proofs of the moments, as well as estimates on 
their rate of growth. 

5.3. General Even Moments. We extend our analysis to consider the general 
even moments. We expect Lemma 15.11 to hold for all general even moments as 
in the case of single palindromic Toeplitz matrices; in other words, we expect in 
the limit as — )■ cxd each matching configuration to contribute equally. Given a 
general even 2rn"^ moment, for each configuration of this moment we have a sys- 
tem of m equations with 2m unknown indices. In the case of single palindromic 
Toeplitz matrices, it is known that we can always choose m + 1 free indices and 
the C-values such that all m equations are satisfied. Furthermore, because C can 
only be {0, ±(A^ — 1)} in the case of single palindromic Toeplitz matrices (all other 
choices give a contribution of size there is only one valid choice for each of 

the remaining m — 1 indices and therefore each configuration contributes 1 to the 
2^th jxioment. However, in the case of highly palindromic Toeplitz matrices there 
are more choices for the m — 1 indices. As the result, the lower order term, which 
was previously negligible, starts contributing to the moment. Nonetheless we con- 
jecture that the contribution from this new term is the same for all configurations 
and sketch a detailed analysis in Appendix |Bl If we could extend Lemma 15.11 to 
the general even 2m*'^ moment, then we would have greatly reduced the complex- 
ity of our moment problem as we would only need to calculate the contribution 
of the completely adjacent matching, and immediately get the same contribution 
from the other (2m — 1)!! — 1 configurations. 

The arguments above, as well as those in the appendices, lead us to make the 
following conjecture. 

Conjecture 5.9 (Configurations Conjecture). Let n be a fixed positive integer, 
and for each N a multiple of 2" consider the ensemble of real symmetric N x N 
palindromic Toeplitz matrices whose first row is 2" copies of a fixed palindrome, 
where the independent entries are independent, identically distributed random vari- 
ables arising from a probability distribution p with mean 0, variance 1 and finite 
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higher moments. Then as N oo each of the (2m — 1)!! configurations for the 
2m^^ moment contribute the same main term. In particular, this means the gen- 
eral 2m}^ moment is just (2m — 1)!! times the contribution from the configuration 
of all adjacent matchings. 

We provide some numerical support for this conjecture in Appendix |X] and some 
theoretical evidence in Appendix [Bl 

5.4. Moment bounds and fat tails. We now extend Theorem 15.71 to matrices 
with greater palindromicity. In doing so, we miss many of the C-vectors that con- 
tribute to these moments, but exact calculations for even a quadruply palindromic 
matrix have proven difficult. The goal is to obtain good enough bounds on the 
moments to deduce properties of the limiting spectral measures. We begin with 
the following lemma. 

Lemma 5.10. Fix n and consider M2m,n, the limit as N —> oo of the average of 
the 2m*'' moment of our ensemble (the set of N x N real symmetric palindromic 
Toeplitz matrices where the first row contains n palindromes) . Then for m suf- 
ficiently large, M2m.,n > (2m — 1)!!, or in other words M2m,n is larger than the 
corresponding moment of the standard normal. If we assume Conjecture \5.9[ then 
we may improve the lower bound to 




M2m,n > 1^-2 ■ (2" - 1) + 2-™" I > ^ c"^ I I . (2m - 1)!!. (5.26) 

As m ^ oo, we have 

M2m,n > (2m -1)!!. (5.27) 

m 

Proof. Let Cc = ^ for b E {1, . . . , 2" — 1}. We determine the contribution to 
the average 2m*'^ moment as — )■ oo from one of the two adjacent matchings. 
That is, consider the core (i.e., the non-zero part) of the corresponding C-vectors 
is (±Cc, tCc, ■ ■ ■ , ±Cc, tCc) and its complement (±(A^ - 1 - Cc),t{N - 1 - 
CJ, . . . , ±{N - 1 - CJ, ^(A^ - 1 - CJ). They contribute 

-2+ (2-;^) + ^) . (5.28) 



2"/ v2" 

The proof of this claim goes back to the observation in Figure [2] that for Cc= 
the number of free / values is A^ — |^ -|- 0(1), whereas if Cc = |^ — 1, then there 
are |^ + 0(1) good / values. Thus the complementary Cc will give the same 
restrictions on the number of / values. 

Moreover, the restrictions on ii from ^ and A^ — 1 — |^ sum to 1. Thus, as 
there are the two cases (plus first or minus first) for each, when we sum them and 
extend the sums back to 0, we have 

--E(';)(^r^-(^)"^(l.)^ (-) 

j even ^ / ^ / ^ / 
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In order to get our lower bound for M2m,n, we repeat this for every value of 
c G {1, 2, . . . , 2"- — 1}. Adding in the contribution from the zero vector, we obtain 



- 2 ■ (2" - 1) + 2-'™ XI ^"^ ' ^^-2°) 



c=l 



which is easily summed for any value of k. 

If we assume Conjecture 15.91 each of the (2m— 1)!! matchings contribute equally, 
and hence 



^2-2"-l 



M,m,n > -2 • (2" - 1) + 2-"^^ ^ c'" ■ (2m - 1)!!. (5.31) 



c=l 



The behavior for large m follows by approximating the sum with an integral. 

If we do not assume Conjecture 15. 9[ we instead use the fact that each config- 
uration contributes at least 1 (this follows from the analysis of the analysis of 
the single palindromic case of |MMSj ) and the analysis above for the case of all 
adjacent matchings. This yields 

M2m,n > + (2m-l)!!-l, (5.32) 

m 

which exceeds (2m — 1)!! for m sufficiently large. □ 

Remark 5.11. We can slightly improve the lower hound in the case when we do 
not assume Conjecture 15.^ hy noting that there are two cases where we have all 
adjacent matchings. Of course, this improvement is negligible in the limit, and it 
is only under the assumption of Conjecture 15.^ that we can significantly improve 
the lower hound. 

We can now turn to an analysis of the properties of the limiting spectral mea- 
sures. Note n = corresponds to the real symmetric palindromic Toeplitz matrices 
studied in |BDJt IMMSj , and n = 1 corresponds to the doubly palindromic Toeplitz 
matrices studied in |MMS] . We now prove Theorem II. 6[ which we restate below 
for the reader's convenience. 



Theorem 11.61 (Fat Tails). Consider the ensemble from Theorem \1.5[ For any 
fixed n > 2, the moments grow faster than the corresponding moments of the 
standard normal. If we additionally assume Conjecture \5.9i then letting M2m,n 
denote the 2m*'* moment of the limiting spectral measure of our ensemble for a 
given n, we have 

M2m,n > (2m -1)!!. (5.33) 

m 

The limiting spectral measure thus has unbounded support, and fatter tails than the 
standard normal ( or in fact any of the known limiting spectral measures arising 
from an ensemble where the independent entries are chosen from a density whose 
moment generating function converges in a neighborhood of the origin). 
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Proof. From Lemma [5.101 we know M2m,n > (2m — 1)!! for m large, and under the 
assumption that Conjecture 15.91 holds we additionally have 

M2m,n > (2m -1)!!. (5.34) 

m 

As n > 1, for m large this is greater than the m^^ moment of the standard normal, 
which is (2m — 1)!!. Thus our limiting spectral measure has unbounded support, 
and more mass in the tails than the standard normal, or in fact, any normal if 
n > 2. To see the last claim, note that if X ~ A^(0, a^) then the 2m^^ moment of 
X is (T^™ ■ (2m — 1)!!, and thus when n > 2 eventually the moment of our ensemble 
is greater than the moment of this normal. □ 

Appendix A. Numerical Methods 

While they can never be accepted as proo^, numerical simulations did much to 
guide our efforts in attacking this problem, and we would not have been successful 
without it, as our naive adaptations of previous works on this subject failed to 
give even remotely accurate predictions. Therefore we give a brief outline of our 
use of these simulations below. 

Initially we primarily used a basic, direct method to approximate the moments 
of the eigenvalue distribution. We first set up a matrix with 2" palindromes 
and chose N so that the matrix had the desired form (every element appears 
exactly 2""*"^ times in the first row). For each moment we used the eigenvalue trace 
lemma to calculate the moment of the eigenvalue distribution for this particular 
matrix, then we averaged over a large number of such random matrices to get an 
approximation for that moment averaged over the ensemble of Toeplitz matrices 
with 2^ palindromes. To get increased accuracy, we simply increased N. 

Our calculations were successful in verifying our predictions for the higher mo- 
ments of the doubly palindromic Toeplitz matrix and for the fourth moment of the 
64-palindrome Toeplitz matrix, supporting our conjectures. Specifically, whenever 
one has involved combinatorial arguments such as the ones above, it is worthwhile 
to numerically test the theory. In Table 1 we present the data from simulating 1000 
real symmetric doubly palindromic 2048 x 2048 Toeplitz matrices, and compare 
the even moments to our predicted values (as expected, the odd moments were 
small). Unfortunately the rate of convergence is slow in N due to the presence of 
large big-Oh constants. 

It is worth noting how slow the convergence is. For example, when we considered 
1000 real symmetric doubly palindromic 96 x 96 Toeplitz matrices, the observed 
second moment was 0.990765, the fourth moment was 4.75209, the sixth was 
45.7965 (for a ratio of 1.22) and the eighth was 737.71 (for a ratio of 1.70). 

While this method was quite useful and accurate for lower moments or for 
a small number of palindromes, for larger values of these quantities the big-Oh 
constants grew quite large, making it computationally prohibitive to simulate a 
representative sample of sufficiently large matrices, and thus leaving us with a 

''When investigating the case of real symmetric pahndromic Toephtz matrices years ago, both 
the authors of |BDJ) and |HM| looked at the numerical simulations and thought the results looked 
Gaussian; it was only after a more careful analysis that the small deviations were isolated. 
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Table 1. Conjectured and observed moments for 1000 real sym- 
metric doubly palindromic 2048 x 2048 Toeplitz matrices. The 
conjectured values come from assuming Conjecture I5.9[ 



Moment 


Conjectured 


Observed 


Observed/Predicted 


2 


1.000 


1.001 


1.001 


4 


4.500 


4.521 


1.005 


6 


37.500 


37.887 


1.010 


8 


433.125 


468.53 


1.082 


10 


6260.63 


107717.3 


17.206 



rather poor estimate of the moments and providing no guide to whether or not 
our formulas were accurate. 

To avoid simulating ever-larger matrices, we instead realized that the average 
2m}^ moment of by matrices in our ensembles should satisfy 

M2^,.^^ = M,^,^ + ^ + ^ + . . . + ^. (A.l) 

Thus, rather than simulating prohibitively large matrices, we could instead simu- 
late large numbers (to increase the likelihood of a representative sample) of several 
sizes of smaller matrices and then perform a least squares analysis to estimate the 
value of M2m,n- We present the data in Table 2. 

In performing the curve fitting, we frequently found big-Oh constants so large 
that it would have been impossible to sample a sufficiently large sample of matrices 
to get an accuracy of within a few percent for the moments. For example, for the 
fourth moment of a Toeplitz matrix with 64 palindromes we found the big-Oh 
constant to be above 30, 000, implying that averages of quite large matrices would 
give an approximation for the fourth moment that would be off by 10, compared 
to a true value of about 128. 

For the doubly palindromic Toeplitz matrices in Table 2, we see that our best 
fit constants for M2m,2 with 2m < 10 are quite close to the values predicted by 
Conjecture 15.91 providing strong evidence supporting it. 

Appendix B. Conjectures for General Configurations 

In Conjecture 15.91 we hypothesize that the main term of the contribution of any 
configuration is equal in the limit. We provide some arguments in support of this. 

To extend Lemma EH] to the general even moment, we introduce some notation 
for a "lift map" , which is a way of relating one configuration of an even moment 
(say one of the (2m — 1)!! configurations of the 2m*'^ moment) to one configuration 
of the next higher even moment (to one of the (2m -|- 1)!! configurations of the 
(2m -|- 2)"'^ moment. If we add a pair of entries to a configuration, this moves 
us from our initial configuration to some configuration of the next even moment. 
There are only two ways to add these entries: adding a pair of adjacent entries or 
adding a pair of non-adjacent entries. 

Conjecture B.l (Configuration Lifting - Adjacent Case). Consider a configura- 
tion of the 2k^^ moment. All configurations of the {2k + 2)"'^ moment obtained 
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Table 2. Observed moments for the ensembles of N x N doubly 
palindromic (i.e., n = 2) Toeplitz matrices. The conjectured values 
are obtained by assuming Conjecture 15. 9[ which means taking the 
contribution from the adjacent matching case and multiplying by 
(2m — 1)!! for the 2m^^ moment. 



N 


sims 


2nd 


4th 


6th 


8th 


10th 


8 


1,000,000 


1.000 


8.583 


150.246 


3984.36 


141270.00 


12 


1,000,000 


1.000 


7.178 


110.847 


2709.61 


90816.60 


16 


1,000,000 


1.001 


6.529 


93.311 


2195.78 


73780.00 


20 


1,000,000 


1.001 


6.090 


80.892 


1790.39 


57062.50 


24 


1,000,000 


1.000 


5.818 


73.741 


1577.42 


49221.50 


28 


1,000,000 


1.000 


5.621 


68.040 


1396.50 


42619.90 


64 


250,000 


1.001 


4.992 


50.719 


858.58 


22012.90 


68 


250,000 


1.000 


4.955 


49.813 


831.66 


20949.60 


72 


250,000 


1.000 


4.933 


49.168 


811.50 


20221.20 


76 


250,000 


1.000 


4.903 


48.474 


794.10 


19924.10 


80 


250,000 


1.000 


4.888 


47.951 


773.31 


18817.00 


84 


250,000 


1.001 


4.876 


47.615 


764.84 


18548.00 


128 


125,000 


1.000 


4.745 


44.155 


659.00 


14570.60 


132 


125,000 


1.000 


4.739 


43.901 


651.18 


14325.30 


136 


125,000 


0.999 


4.718 


43.456 


637.70 


13788.10 


140 


125,000 


1.000 


4.718 


43.320 


638.74 


14440.40 


144 


125,000 


1.001 


4.727 


43.674 


647.05 


14221.80 


148 


125,000 


1.000 


4.716 


43.172 


628.02 


13648.10 


Conjectured 




1.000 


4.500 


37.500 


433.125 


6260.63 


Best Fit M2m,2 




1.000 


4.496 


38.186 


490.334 


6120.94 



by adding a pair of adjacent entries to this configuration contribute equally to the 
(2k + 2)'"''^ moment. 

We comment on the above conjecture; see Figure H] for an example. Let 

(. . . , Opq, . . . , Ojj, Ojfc; • • • (^Imi • • • ) 

be a tuple from one of the (2A; — 1)!! configurations of the 2k^^ moment; for brevity's 
sake we call this configuration (1). We let 

(. . . , Clpg, . . . , fljj, Ctjo) O'osi (^ski • • • ; (^Imi • • • ) 

be the new tuple obtained by adding the pair of adjacent, matched entries a jo = 
Qos] we denote this by configuration (2). Let fl2k be the set of all tuples that work 
for configuration (1) and ^l2k+2 be the set of all tuples that work for configuration 
(2). We define a "lift map" F : fiafc -> ^2k+2 by 

F ((. . . , Ctpg, . . . , fljj, Ojr'fe; • • • ; (^Imi • • • )) (• • • 5 ^pqi ■ ■ ■ i (^iji (^joi ^osi (^ski ■ ■ ■ i ^Imy ■ ■ ■ j 

(B.l) 

Note F maps each index in (. . . , a^g, . . . , ajj, a^-fc, . . . , aim-, . . . ) to itself and inserts 
a new index s = j — B + B where B is the value of C corresponding to the pair 
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Figure 4. Moment Lifting by adding a pair of adjacent entries 



of entries {ajk — Op^), and B' is any value of C such that s E {1, . . . , N} and 
{B — B') is a, vahd value of C. The system of equations corresponding to the 
tuples (. . . , Qpq, . . . , Qij, Qjk, ■ ■ ■ , aim, ■ ■ ■ ) IS givcu as follows: 



I 



l — m — —{i—j)-\-A 

P-q = -U-k) + B 



(B.2) 



l — m — 




p-q = 


-is 


j-o = 


-io 



Under the map F, we obtain a new tuple (. . . , Upg, . . . ,aij, ajo, aos, cisk, ■ ■ ■ , dim, ■ ■ ■ ) 
satisfying the system of equations 



(B.3) 



Except for the two equations p—g = —{s — k)+B' a,ndj — o— —{o~s) + {B — B'), 
every other equation of configuration (1) is preserved under F and therefore still 
holds in configuration (2). Furthermore, since both B and {B - B) are valid 
choices of the C value by the construction of F, the two equations p — q — —{s — 
k) + B and j — o = —{o — s) + {B — B ) are also valid. Thus 



(. . . , a 



pq, 



, dsk, ■ ■ ■ 1 dim, ■ ■ ■ ) ^ ^2fc+2- 



(B.4) 
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Based on the analysis above, the mapping F only depends on the choice of one 
index j and one C value B from the original configuration (1). First, we can take 
B to be any possible value of our C values since in order to obtain all the tuples 
of configuration (1), we need to sum over all possible combinations of the C's that 
work for the system of equations corresponding to configuration (1). 

Second, for any configuration in the 2k^^ moment, we have 2k indices (unknown 
variables) and k equations with the last equation linearly dependent on the rest. 
Thus we must have at least two completely free indices that can take on any value 
between 1 and N. We specify the two completely free indices by choosing the very 
first entry of the tuple, which obviously can be any vertex of the configuration. 
Starting with any vertex, if there are some m satisfying tuples with aij being our 
first entry, then there must be exactly m tuples with aji being our first entry 
since the matrix is symmetric over the main diagonal. Therefore, the number 
of resulting tuples from the map F is the same regardless of where we add the 
adjacent pair. 

The biggest obstacle left is whether the map F can reach every possible tuple 
of configuration (2). Given a tuple (. . . , apg, . . . , a^j, ajo, aos, ask, ■ ■ ■ , o-im, • • • ) of 
the configuration (2), the inverse map F"^ simply substitutes the equation arise 
from the added pair of adjacent entries {ajo = jos) j = s + into the equation 
arise from the entries {ask = jpq) p — q = —{s — k) + Ci. In order for this tuple 
to be reachable by the map F, it suffices that C2 + C3 is a valid C value. It 
is not a problem for the 4^^ and 6^^ moment since the sum of all the C values 
corresponding to any tuple must equal 0. However, at higher even moment, there 
is possibility that C2 + C3 is not a valid C value. So there might be tuples of the 
configuration (2) that the map F can not reach. If those tuples only contribute to 
some lower order term then all configurations of the {2k + 2)'^'^ moment obtained 
by adding a pair of adjacent entries to a given configuration of the (2^)**^ moment 
would contribute equally. 

Conjecture B.2 (Configuration Lifting - Nonadjacent Case). Consider a config- 
uration of matchings for the 2k^'^ moment. All configurations at the 2k^^ moment 
obtained by adding a pair of non- adjacent entries contribute equally to the (2/c+2)"'' 
moment. 

We provide some arguments in support of our claim; see Figure |5]for an example. 
Let (. . . , aij, ttjk, . . . , aip, apg, . . . ) be a tuple of a configuration for the 2/c*^ moment; 
denote this by configuration (1). We let (. . . , a^j, ajo, aok, ■ ■ ■ , aip, aps, a^g, . . . ) be 
the new tuple obtained by adding the pair of entries As before, let 

^2?: be the set of all tuples that work for configuration (1) and ^'2k+2 be the 
set of all tuples that work for the configuration (2). We define the "lift map" 
FF : figfc ~^ ^2k+2 this case by 

F {{. . . , fljj, ajk, . . . , aip, apq, . . . )) (. . . , a^j, ajo, aok, ■ ■ ■ aip, aps, agq, • • • ) (B.5) 

such that F' maps every index in (. . . , a^j, ajm, • • • , o/p, apg, . . . ) to itself and adds 
two new indices = j + B — B and s = p + D — D where B and D are the C 
values associated with the pairs containing ajm and apg respectively. Also, B and 
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Figure 5. Moment Lifting by adding a pair of non-adjacent entries 

D' are any value of C such that o, s E {1, . . . , A^} and {D' + B' — D — B) is a 
value of C. 

For the tuple (. . . , aij, ajm, • • • , a^p, ctpg, . . .) we have the following system of equa- 
tions: 

'io-jo = -(^ - j) + A 
Jo-ko = -{j-k) + B 
< lo-Po = -{l-p) + C (B.6) 
Po - go = -{p -q) + D 
k ■ ■ ■ 

Under the map F\ we obtain a new tuples (. . . , ajj, a jo, aok, 
satisfying the system of equations: 

'^0 - Jo = -(^ - j) + A 
jo - ko = -(o - k) + B' 
^ k-po = -{l-p) + C 
Po-qo = -{s -q) + D' 
j-o = -{p-s) + {D' + B' -D-B) 

k ' ' ' 

Similar to the analysis in Conjecture IB.lt ^ill equations except for jo ~ ^o = 
-{o-k) + B' andpo-^o = -{s-q) + D' and j-o = -{p- s) + {D' + B' - D - B) 
are preserved under the map F so they still hold. Furthermore, since B , D 



. . . , Clip, Ojps, Clsq, ■ ■ ■ 



(B.7) 
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and {D' + B' ~ D — B) are all valid choices of C, the other three equations also 
hold true. Lastly, the existence of at least two completely free indices allow us 
to choose them to be the first index in ajm and the second index in aip. Thus, 
following the same line of argument in Conjecture IB.lt can always choose the 
two free indices such that the number of tuples resulting from the map F (^2/=) 
are the same regardless of where we add the non-adjacent pair. 

Corollary B.3. Given any configuration, we can replace one of its adjacent pairs 
by another adjacent pair, and similarly for non- adjacent pairs, without changing 
its contribution to the corresponding moment. 

Given any configuration at the {2k + 2)°^^ moment. Corollary IB. 31 allows us to 
repeatedly replace adjacent pairs with other adjacent pairs, and similarly for non- 
adjacent pairs. By iterating this process, we can move any configuration down 
to the completely adjacent configuration (which contains only adjacent pairs of 
entries); see for instance Figure El Given any configuration, whenever there is a 
crossing between two pairs of entries, which implies those entries are non-adjacent, 
we can eliminate the crossing by replacing one of these pairs of entries with another 
pair of non- adjacent entries. 




Figure 6. Eliminating crossings in a given configuration 

Iterating this process, we end up with a configuration that has no crossing be- 
tween any of its pairs of entries. We then keep replacing the pair of adjacent entries 
at the two ends of the configuration with another pair of adjacent entries in the 
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middle of the configuration until we obtain the completely adjacent configuration; 
see Figure [3 




Figure 7. Transformation to completely adjacent configuration 

Based on Corollary IB.3I and Conjecture 15.11 it follows by induction that every 
configuration at any even moment would contribute equally. 
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